Search for: All records

One-dimensional stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleosynthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundary that induce chemical and angular momentum transport, as well as modify the thermal structure of the star. One such process is called convective penetration, wherein vigorous convection extends beyond the nominal convective boundary and alters both the composition and thermal structure. In this work, we incorporate the process of convective penetration in stellar evolution calculations using the stellar evolution software instrumentmesa. We implement convective penetration according to the description presented by Anders et al. to to calculate a grid of models from the pre-main sequence to helium core depletion. The extent of the convective penetration zone is self-consistently calculated at each time step without introducing new free parameters. We find both a substantial penetration zone in all models with a convective core and observable differences to global stellar properties such as the luminosity and radius. We present how the predicted radial extent of the penetration zone scales with the total stellar mass, age, and metallicity of the star. We discuss our results in the context of existing numerical and observational studies.

 

Massive stars die in catastrophic explosions that seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion’s character and the remnant mass depend on models of the star’s evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present three-dimensional simulations of massive star convection extending from the star’s centre to near its surface, with realistic stellar luminosities. Using these simulations, we predict the photometric variability due to convectively driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, and thus the red noise must be generated by an alternative process.

 

We revisit the tidally excited oscillations (TEOs) in the A-type main-sequence eccentric binary KOI-54, the prototype of heartbeat stars. Although the linear tidal response of the star is a series of orbital-harmonic frequencies which are not stellar eigenfrequencies, we show that the non-linearly excited non-orbital-harmonic TEOs are eigenmodes. By carefully choosing the modes which satisfy the mode-coupling selection rules, a period spacing (ΔP) pattern of quadrupole gravity modes (ΔP ≈ 2520–2535 s) can be discerned in the Fourier spectrum, with a detection significance level of $99.9{{\ \rm per\ cent}}$. The inferred period spacing value agrees remarkably well with the theoretical l = 2, m = 0 g modes from a stellar model with the measured mass, radius, and effective temperature. We also find that the two largest-amplitude TEOs at N = 90, 91 harmonics are very close to resonance with l = 2, m = 0 eigenmodes, and likely come from different stars. Previous works on tidal oscillations primarily focus on the modelling of TEO amplitudes and phases, the high sensitivity of TEO amplitude to the frequency detuning (tidal forcing frequency minus the closest stellar eigenfrequency) requires extremely dense grids of stellar models and prevents us from constraining the stellar physical parameters easily. This work, however, opens the window of real tidal asteroseismology by using the eigenfrequencies of the star inferred from the non-linear TEOs and possibly very-close-to-resonance linear TEOs. Our seismic modelling of these identified eigen g-modes shows that the best-matching stellar models have (M ≈ 2.20, 2.35 M⊙) and super-solar metallicity, in good agreement with previous measurements.

 

ABSTRACT The first bright objects to form in the Universe might not have been ‘ordinary’ fusion-powered stars, but ‘dark stars’ (DSs) powered by the annihilation of dark matter (DM) in the form of weakly interacting massive particles (WIMPs). If discovered, DSs can provide a unique laboratory to test DM models. DSs are born with a mass of the order of M⊙ and may grow to a few million solar masses; in this work we investigate the properties of early DSs with masses up to $\sim \! 1000 \, \mathrm{ M}_\odot$, fueled by WIMPS weighing 100 GeV. We improve the previous implementation of the DM energy source into the stellar evolution code mesa. We show that the growth of DSs is not limited by astrophysical effects: DSs up to $\sim \!1000 \, \mathrm{ M}_{\odot }$ exhibit no dynamical instabilities; DSs are not subject to mass-loss driven by super-Eddington winds. We test the assumption of previous work that the injected energy per WIMP annihilation is constant throughout the star; relaxing this assumption does not change the properties of the DSs. Furthermore, we study DS pulsations, for the first time investigating non-adiabatic pulsation modes, using the linear pulsation code gyre. We find that acoustic modes in DSs of masses smaller than $\sim \! 200 \, \mathrm{ M}_\odot$ are excited by the κ − γ and γ mechanism in layers where hydrogen or helium is (partially) ionized. Moreover, we show that the mass-loss rates potentially induced by pulsations are negligible compared to the accretion rates. 

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The newauto_diffmodule implements automatic differentiation inMESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection inMESAwith a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthenMESA’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars inMESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates toMESA’s software infrastructure that enhance source code development and community engagement.

 

ABSTRACT Strongly magnetic B-type stars with moderately rapid rotation form ‘centrifugal magnetospheres’ (CMs) from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A long-standing question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic ‘centrifugal breakout’ (CBO) events, wherein magnetic tension can no longer contain the built-up mass. We argue here that recent empirical results for Balmer-α emission from such B-star CMs strongly favour the CBO mechanism. Most notably, the fact that the onset of such emission depends mainly on the field strength at the Kepler radius, and is largely independent of the stellar luminosity, strongly disfavours any drift/diffusion process, for which the net mass balance would depend on the luminosity-dependent wind feeding rate. In contrast, we show that in a CBO model, the maximum confined mass in the magnetosphere is independent of this wind feeding rate and has a dependence on field strength and Kepler radius that naturally explains the empirical scalings for the onset of H α emission, its associated equivalent width, and even its line profile shapes. However, the general lack of observed Balmer emission in late-B and A-type stars could still be attributed to a residual level of diffusive or drift leakage that does not allow their much weaker winds to fill their CMs to the breakout level needed for such emission; alternatively, this might result from a transition to a metal–ion wind that lacks the requisite hydrogen. 

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729M_⊙carbon–oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of²²Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6−7 significant figures of precision. Usually white dwarf models without²²Ne are fit to the observed periods and other properties. The rms residuals to the ≃150−400 s low-order g-mode periods are typically in the range ofσ_rms≲ 0.3 s, for a fit precision ofσ_rms/P≲ 0.3%. We find average relative period shifts of ΔP/P≃ ±0.5% for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of ΔP/Pdepending on the specific g-mode, abundance of²²Ne, effective temperature, and the mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when²²Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of²²Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including²²Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model-fitting processes.